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Abstract 

Formulas for stable differentiation of piecewise-smooth functions are given. The 
data are noisy values of these functions. The locations of discontinuity points and 
the sizes of the jumps across these points are not assumed known, but found stably 
from the noisy data. 

1 Introduction 

Let / be a piecewise-C^([0, 1]) function, < Xi < X2 < ■ ■ ■ < xj,l < j < J, are 
discontinuity points of /. We do not assume their locations xj and their number J known 
a priori. We assume that the limits f{xj ± 0) exist, and 

sup If^^^^x)] < Mrn, m = 0,1,2. (1.1) 

X^Xj,l<j<J 

Assume that fs is given, ||/ - fs\\ := sup^^^. ^^j^j \ f - fs\ < 5, where fs G L°°(0, 1) are 
the noisy data. 

The problem is: given {fs, 6}, where 6 G (0, 6o) and 6o > is a small number, estimate 
stably f , find the locations of discontinuity points xj of f and their number J, and estimate 
the jumps pj := f{xj + 0) — f{xj — 0) of f across Xj, I < j < J. 

A stable estimate Rsfs of /' is an estimate satisfying the relation lim^^o WRsfs — f'\ \ = 

0. 

There is a large literature on stable differentiation of noisy smooth functions (e.g., see 
references in 0), but the problem stated above was not solved for piecewise-smooth func- 
tions by the method given below. A statistical estimation of the location of discontinuity 
points from noisy discrete data is given in PJ. 
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The following formula was proposed originally (in 1968, see references in j2]) for stable 
estimation of f'{x), assuming / G C^([0, 1]), M2 7^ 0, and given noisy data fs: 

^^^^^^ A(x + Mi)W.(x-Mi)) _,(,)^, (H)', h(S)<.<i-h(S). (1.2) 

and 

WRsfs-fW <^/2M^5:=e{5), (1.3) 

where the norm in (1.3) is L°°(0, 1)— norm. Moreover, (cf i2]), 

inf sup \\Tfs- f\\>e{5), (1.4) 

where T : L°°(0, 1) L°°{0, 1) runs through the set of all bounded operators, K{M2, 5) : = 
{f '■ 1 1 /"I I < M2, \\f — fs\\ < 5}. Therefore estimate ()1.2|1 is the best possible estimate of 
/', given noisy data fs, and assuming / G K{M2, S). 

In j2] this result was generalized to the case / G K{Ma,S), ||/^"^|| < Ma, 1 < a < 2, 
where ||/(")|| := ||/|| + ||/' || +sup^_^, ^^^P^' 1 < « < 2, and f^""^ is the fractional-order 
derivative of /. 

The aim of this paper is to extend the above results to the case of piecewise-smooth 
functions. In Section 2 the results are formulated, and proofs are given. In Section 3 the 
case of continuous piecewise-smooth functions is treated. 



2 Formulation of the result 

Theorem 1. Formula ()1.2|) gives stable estimate of f on the set Sg '■= [h{5),l — 
^('^)] \ Uj=i(^j ~ h(5),Xj + h(6)), and ()1.3j) holds with the norm \\ ■ \\ taken on the 

set Ss- Assuming M2 > and computing the quantities fj := where 
1 

h := h{6) := (^'^^^ > 1 ^ i < for sufficiently small S, one finds the location of 

discontinuity points of f with accuracy 2h, and their number J. Here [^] is the integer 
smaller than ^ and closest to ^. The discontinuity points of f are located on the intervals 
{jh — h,jh + h) such that \fj\ ^1 for sufficiently small 6, where e{6) is defined in (1.3). 
The size pj of the jump of f across the discontinuity point Xj is estimated by the formula 
Pj ~ fsiJh + h) — fs{jh — h), and the error of this estimate is 0{y/6). 

Let us assume that miUjPj : = p ^ h{d), where ^ means "much greater than". Then 
Xj is located on the j-th interval [jh — h,jh + h], h := h{6), such that 

f5Uh + h)~fs{jh-h) 



I/. 



2h 



» 1, (2.1) 



so that Xj is localized with the accuracy 2h{S). More precisely, \fj\ > ^^^^'^^'^^ ^^^^ — ^ 



2h /i' 



and ^ = 0.5e{6), where e{5) is defined in (1.3). One has 
\f{jh + h)-fUh-h)\>p,-\f{jh + h)-f{x, + 0)\-\fUh-h)-f{x,-0)\>p,-2M,h. 
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Thus, 



|/,| > H - Ml - 0.5e{S) = ci^ - C2 » 1, 



where Ci := and Ca := Mi + 0.5£(5). 

The jump pj is estimated by the formula: 

Pj^[fsUh + h)-fsUh-h)], (2.2) 

and the error estimate of this formula can be given: 



\Pj - [fsijh + h)- fsijh -h)]\<26 + 2M,h = 26 + SMi^ ^ = 0{V6). (2.3) 

Thus, the error of the calculation of pj by the formula pj ~ /^(j/i + h) — fs{jh — h) is 
0{6^) as 6^0. 

Proof of Theorem 1. If x G Ss, then using Taylor's formula one gets: 

\{Rsfs){x)-f{x)\<^- + ^. (2.4) 

Here we assume that M2 > and the interval {x — h{6), x + h{6)) C Ss, i.e., this interval 
does not contain discontinuity points of /. If, for all sufficiently small h, not necessarily 
for h = h{S), inequality (2.4) fails, i.e., if |(i?5A)(x) - f'{x)\ > J + ^ for all sufficiently 
small h > 0, then the interval {x — h, x + h) contains a point Xj ^ Ss, i.e., a point of 
discontinuity of / or /'. This observation can be used for locating the position of an 
isolated discontinuity point xj of / with any desired accuracy provided that the size \pj\ 
of the jump of / across xj is greater than 46, \pj\ > 46, and that h can be taken as small 
as desirable. Indeed, if Xj e {x — h, x + h), then we have 26 + Ma/i^ < \ fs{x + h) — fs{x — 
h)—2hf'{x)\ = \pj + 0{h'^)±26\. The above estimate follows from the relation \ fs{x + h) — 
fs{x-h)-2hf'ix)\ = \fix+h)-fix,+0)+Pj+fix,-0)-fix-h)±26\ = \p,+0{h^)±26\. 
Here |p± &|, where 6 > 0, denotes a quantity such that \p\ — b< \p±b\ < \p\ + b. Thus, if 
h is sufficiently small and \pj \ > 46, then the inequality 26 + Ma/i^ < \pj + 0{h^) ± 26\ can 
be checked, and therefore the inclusion h,x + h) can be checked. Since h > is 

arbitrarily small in this argument, it follows that the location of the discontinuity point 
Xj of / is established with arbitrary accuracy. Additional discussion of the case when a 
discontinuity point Xj belongs to the interval {x — h{6), x + h{6)) will be given below. 

Minimizing the right-hand side of (|2.4p with respect to h yields formula (|1.2p for the 
minimizer h = h{6) defined in (1.2), and estimate (jl.3|) for the minimum of the right-hand 
side of (El- 

If p ^ h{6) , and ()2.1|) holds, then the discontinuity points are located with the accu- 
racy 2h{6), as we prove now. 

Consider the case when a discontinuity point Xj of / belongs to the interval {jh — 
h,jh + h), where h = h{6). Then estimate ()2.2|) can be obtained as follows. For jh — h < 
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< jh + h, one has 

\f{xj + 0) - fix, - 0) - fsUh + h) + fsijh -h)\< 2S+ 
+ \f(x, + 0)-fUh + h)\ 

+ \f{xj - 0) - f{jh -h)\<26 + Mi2h, h = h{6). 

This yields formulas (2.2) and (2.3). Computing the quantities fj for 1 < j < [^], and 
finding the intervals on which (2.1) holds for sufficiently small d, one finds the location 
of discontinuity points of / with accuracy 2h, and the number J of these points. For 
a small fixed 6 > the above method allows one to recover the discontinuity points of 
/ at which \fj\ > ^ - f - Mi > 1- This is the inequality (2.1). If h = /i(5), then 
I = 0.5e{5) = 0{V6), and \2hfj-pj\ = 0{VS) as 5 provided that M2 > 0. Theorem 
1 is proved. □ 
Remark 1: Similar results can be derived if ||/*^°^||icx)(s^) := < Ma, 1 < a < 2. 

In this case h = h{5) = CaS'^, where Cq = Rsfs is defined in (jl.2|) . and the 

error of the estimate is: 

a-l 

\\R5fs-f'\\ss<aMa^ (^] " 6 



a-l 



The proof is similar to that given in Section 3. It is proved in [2J that for C"-functions 
given with noise it is possible to construct stable differentiation formulas if a > 1 and it 
is impossible to construct such formulas if a < 1. The obtained formulas are useful in 
applications. One can also use L^-norm on Ss in the estimate Wf^^^Ss < Ma (cf. 0). 

Remark 2: The case when M2 = requires a special discussion. In this case the last 
term on the right-hand side of formula (2.4) vanishes and the minimization with respect 
to h becomes void: it requires that h be as large as possible, but one cannot take h 
arbitrarily large because estimate (2.4) is valid only on the interval {x — h,x + h) which 
does not contain discontinuity points of /, and these points are unknown. If M2 = 0, then 
/ is a piecewise-linear function. The discontinuity points of a piecewise-linear function 
can be found if the sizes \pj\ of the jumps of / across these points satisfy the inequality 
\pj\ >> 26 + AMih for some choice of h. For instance, iih = -^^^ then 25+AMih = 45. So, 
if \pj\ » 46, then the location of discontinuity points of / can be found in the case when 
M2 = 0. These points are located on the intervals for which \ fs{jh+h) — fs{jh—h)\ » 46, 
where h=^. 

The size \pj\ of the jump of / across a discontinuity point Xj can be estimated by 
formula (2.2) with h = and one assumes that Xj G {jh — h,jh + h) is the only 

discontinuity point on this interval. The error of the formula (2.2) is estimated as in the 
proof of Theorem 1. This error is not more than 26 + 4Mih = 46 for the above choice of 

One can estimate the derivative of / at the point of smoothness of / assuming M2 = 
provided that this derivative is not too small. If M2 = 0, then / = ajx + bj on every 
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interval between the discontinuity points xj, where aj and bj are some constants. If 

{jh-h,jh+h) C Aj, and fj := M^'^+^'^-Mjh-h) ^ ^^^^ \fj-aj\ < f Choose = ^, where 
t > is a parameter, and Mi = maxj \ aj\. Then the relative error of the approximate 
formula aj ~ /-,• for the derivative /' = aj on equals to ^•^^^ < Thus, if, e.g., 

> ^ and t = 10, then the relative error of the above approximate formula is not 
more than 0.1. 



3 Continuous piecewise-smooth functions 

Suppose now that ^ e {mh — h, mh + h) , where m > is an integer, and ^ is a point at 

which f is continuous but f'{C) does not exist. Thus, the jump of / across ^ is zero, but 
^ is not a point of smoothness of /. How does one locate the point ^ ? 

The algorithm we propose consists of the following. We assume that M2 > on Ss- 
Calculate the numbers fj := Mf'+''^~Mjh-h) |^^^^ _ _ 2, . . . , /i = h{S) = 

Y^^. Inequality (1.3) implies fj — £{S) < f'{jh) < fj + s{S), where £{S) is defined in 
(1.3). 

Therefore, if \fj\ > £{5), then sign fj — sign f'{jh). 
One has: 

7-f<i/,«-/,i<^+f, 

where f = 0.5^(5) and J := | fUh+2h)-f{jh)-f{jh+h)+f{jh-h) | _ -y^.^g Taylor's formula, one 
derives the estimate: 

0.5[Ji - e{5)] <J< 0.5[Ji + e{5)], (3.1) 

where Ji := \f'{jh + h) - f'{jh)\. 

If the interval {jh-h,jh + 2h) belongs to Ss, then Ji = \f{jh + h) - f'{jh)\ < Ma/i = 
s{d). In this case J < e{5), so 

\fj+i-fj\<le{6) if {jh-h,jh + 2h)GSs. (3.2) 

Conclusion: if \ fj+i — fj\ > then the interval {jh — h,jh-\- 2h) does not belong 

to Ss, that is, there is a point ^ G {jh — h,jh + 2h) at which the function f is not twice 
continuously differentiable with |/"| < M2. Since we assume that either at a point ^ 
the function is twice differentiable, or at this point f does not exist, it follows that if 
\fj+\ — fj\ > j^iS), then there is a point ^ e {jh — h,jh + 2h) at which f does not exist. 

If 

f if 3+1 < 0, (3.3) 

and 

min(|/,+i|,|/,|)>£(<5), (3.4) 

then (3.3) implies f'{jh)f'{jh + h) < 0, so the interval {jh,jh + h) contains a critical 
point of /, or a point ^ at which /' does not exist. To determine which one of these 
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two cases holds, let us use the right inequality (3.1). If ^ is a critical point of / and 
^ ^ {jh,jh + h) C Ss, then Ji < e{6), and in this case the right inequality (3.1) yields 

<e(5). (3.5) 

Conclusion: If (3.3)-(3.5) hold, then ^ is a critical point. If (3.3) and (3.4) hold and 
l/j+i — fj\ > s{6) then ^ is a point of discontinuity of f. 

If ^ is a point of discontinuity of /', we would like to estimate the jump 

i^:=i/'(e+o)-/'(e-o)|. 

Using Taylor's formula one gets 

.f,+i-fj = ^ ±3.56(5). (3.6) 

The expression A = B ± b, b > 0, means that B — b<A<B + b. Therefore, 

P = 2(/,+i-/,)±75(<5). (3.7) 
We have proved the following theorem: 

Theorem 2. If G (jh — h,jh + 2h) is a point of continuity of f and — fj\ > ^^{5), 
then ^ is a point of discontinuity of f . If ()3.3|) and fl3.4p hold, and \ fj+i — fj\ < s{5), 
then ^ is a critical point of f . If fl3.3|) and ()3.4p hold and \fj+i ~ fj\ > ^{5), then 
^ G {jh,jh + h) is a point of discontinuity of f. The jump P of f across ^ is estimated 
by formula (3.7). 
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